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We calculate the electron-mediated exchange interaction between two domain walls in magnetic wires. This corre- 
sponds to the equilibrium regime and, therefore, the interaction can be additionally controlled by an electric current. 
The exchange interaction is long ranged and oscillates as a function of the distance between the walls. It also depends 
oscillatory on the polarization angle of the walls, having the maximum value for collinear polarization. 
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1 Introduction Magnetic domain walls (DWs) have 
attracted a lot of attention as important elements of new 
magnetoelectronic devices OEEI. In particular, it was 
demonstrated recently that they can be used in a new type 
of memory device (the racetrack memory) effectively con- 
trolled by an electric current [4 5|. On the other hand, the 
DW can be viewed as a kind of local imperfection in an 
ordered magnetic system, like an impurity or defect. The 
substantial difference is that DWs can move and therefore 
they can be relatively easily put into motion by an external 
field or electric current, and also by an interaction between 
them. 

A great amount of theoretical and experimental work 
was dedicated to studying the resistance of DWs, the 
current-induced spin torque, the dynamics of DW motion, 
and other effects related to a single DW strongly coupled 
to the electron system {SEE). When the density of DWs 
increases, it is important to include into these considera- 
tions the effects of their interaction. It has already been 
demonstrated |8 ') 10 1 Q that an electric current in a mag- 
netic wire with DWs influences the DW interaction, so that 
by using a current the DW coupling can be controlled, and 
the dynamics of strongly coupled DWs can be affected. 

In this work we consider in detail the indirect exchange 
coupling between the DWs in equilibrium. Essentially, the 
analysis of such interaction is the first necessary step to 
understand the basic mechanisms of DW interactions. 
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2 Model Let us consider the following Hamiltonian, 
describing an electron gas coupled by exchange energy M 
to a textured magnetization which changes its orientation 
at each DW: 

H = — -M<r. n (r), (1) 

2m 

where the unit vector n(r) determines the magnetization 
orientation, m is the carrier effective mass, and we take 
units with h = 1, We assume the ferromagnetic wire (or 
ribbon) to be orientated along the x axis and consider two 
DWs, labeled as 1 and 2. For definiteness, we also assume 
that at x — > — oo and x — > +oo the magnetization is along 
unit vector no = (0, 0, 1), and within each of the DWs the 
vector n is rotated by an angle tt around the unit vector tj, 
where the index % = 1, 2 refers to DW 1 or 2. We focus here 
on the case of transverse DWs, the results for a vortex DW 
should look qualitatively similar. Thus, the dependence of 
n(x) within DWs 1 and 2 is described as 



n 1 (x) = e^ ltl - L n , n 2 (x) = _ e «^t a -L ^ (2) 

where L is the matrix of moment L — 1 and tpi (x) changes 
from to tt when x crosses the i-th DW. Correspond- 
ingly, the transformation of the spinor wave function ip to 
the local frame with the magnetization along the axis z is 



ip — > Tip = e 2 



ip. The location Xi and the vector t; 



determine the DW state. Note that for the electron motion 
the DW can be always considered as static even for DWs 
moving along the wire. 
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We use this transformation to the local frame 11211131 
FBI , in which the vector of each DW is oriented in the 
same direction along the global axis z. After this transfor- 
mation the electron gas is in the homogeneous magnetiza- 
tion M but there appears the gauge potential related to the 
local transformation, A;(.t) = i((p' i /2) (tj • <r). 

We assume the DW width A to be much larger than the 
electron wavelength at the Fermi surface, kp\^> 1, which 
is the typical condition for metallic ferromagnets. Then the 
transformed Hamiltonian is (summing over i) 



H = y 



fc 



2m 



2m 



- Ma, - 



Pi 



h.c. 



, (3) 



where Pi(x) = (p'(x — Xi)/2m, p x = —id/dx and X{ is 
the point where the i-th DW is located. 

The exchange interaction energy can be found using a 
RKKY approach with the DW-induced perturbation local- 
ized in the vicinity of the points x\ = and X2 = R. Using 
Eq. (3) we find 



E; 



a cs Rc Tr 



d 2 k dp 

' ^ dx' dx" pt(x') fo(x") 



x(t 



1-<T 



(2tt) 2 2tt 
dG ke (x' -x") 



dx' dx" 
where a cs is the cross-section of the DW and 

rcc dk 



G kE (x) = 



2tt 



diag (Gfcet ) G kel) 



(5) 



is the Green function of electrons in a homogeneous mag- 
netization field: G^ a = e — Eka — k^/2m + fi + i5 sgn(e), 
where \i is the chemical potential. We also denoted 
£fc = (ky + fc 2 )/2m and £kf,i = £ k T M. Now we take 
k a = +y/2m(e — e^a + A f ) + ^ sgne. Then defining 
Cfeo- = e k<7 — M we have the conditions that if e — £ka > 
and e < then k a lies in the lower half-plane, otherwise 
it lies in the upper half plane. This allows us to directly 
calculate Gk £ (x). 

Calculating the Green functions and their derivatives 
in coordinate representation and substituting them into 
Eqn. (H]i we obtain 



Ej 



-ff r ,m Re Tr 



d 2 k de 



dx dx" 



(2tt) 2 2tt 
x ^(x') I3 2 (x" ) ■ a) P(t 2 - a) P, 



(6) 



where we denoted 

P = diag(./- fet e 4fc ^"-*') + (i _ /fet ) e ifet(x"-x') S gn(s) ) 
f H e- ik ^ x '- x ") + (1 - / H ) e -<Mx'-x")«n(*)) (?) 

and ff-a = f{£ka) is the Fermi-Dirac function at T — > 0. 
The traces over the matrices can be performed immedi- 
ately. Assuming that the first DW points into the y-axis, 
then we can take ti = x and t 2 = x cos 9 + y sin 8, so that 




Figure 1 The interaction energy as a function of inter DW 
distance R, for A = lOA^ and two different angles be- 
tween the DWs. See text for details. 



9 is the angle between the two domain wall polarizations. 
At T = this yields 

Tr (t x • u) P (t 2 • a) P = 2 cos(0) 

x [/ fe -|-e _ ^ x ' _x "^ fet+fe ' 1 -^ + e -j(x'-x")(A;t+fcj.sgnE) 
+ (1 - fkl)e~ l{x '-x")(feT+fc4.)sgn £ j _ (g) 

For the /3 function we use j3i (x) = sech(a; / L) /2mL, then 
we can calculate the integrals over dx' and dx" using 



7i Ek 



dx' f3i(x') e Tikx = — sech 
n ; 2m 2 



(9) 



for i = 1 , 2 and i?i = 0, R 2 — R- Now we can switch the 
fc-integrals to polar coordinates and perform the angular 
integral, and then substituting e' = k 2 /2m and rescaling 

e = e — e' + /i we obtain 



Eint(9, R) — ^rz R 



?,2 



de' 

J M 



dee iRkl 



x scch 



TrLjk! + k 2 ) r iR k 2+e iRk 2 agn(e+e'- f j,) 

2M f&e^sech' -Hh + h) 
Jm 2 



(10) 



where k\,2 = y2m(e ± M). 

The dependence of the exchange interaction on the dis- 
tance for the two polarization angles 9 is shown in Fig. Q] 
Here we use the following parameters: the Fermi wave- 
length A f = 0.367 nm, M = 36 meV, <r cs = 20 x 20 nm 2 , 
and A = IOAf- The magnitude of the interaction depends 
strongly and non-mono tonic ally on the DW width, A, and 
the magnetization strength, M. They are related by the 
strengths of the anisotropy and the exchange energy in the 
material. Halving the width of the DW to A = 5Af, we 
already see a much larger effect, see Fig. [2] 

The obtained results show that the interaction of two 
DWs is long-ranged and strongly oscillating. The particu- 
lar behaviour around R re 4A is caused by the change of 
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Figure 2 The interaction energy as a function of inter DW 
distance R, for 9 = 7r/4 and A = 5Ap\ See text for details. 

sign of the envelope function of the energy. I.e. the enve- 
lope function is also not a monotonic function of R. But 
the most important is that the DW interaction depends on 
both the distance between the walls and the DW polariza- 
tion determined by the angle 9. It means that if we put 
DWs located at a certain distance R, which corresponds to 
the energy-favorable collinear mutual polarization of DWs, 
then there is another location very close in R with the anti- 
collinear orientation of DWs, with almost the same energy. 
Considering the energy of the system as a function of dis- 
tance R, we find the correspondence to a classical "parti- 
cle" in the oscillating potential profile, so that the neighbor- 
ing positions of this particle in the minima of the potential 
describe the up and down states of one of the DWs with 
respect to the other one. As the amplitude of interaction in- 
creases with decreasing R at small distances (see Figs. 1 
and 2), the DWs are effectively attracted to each other. 

One can assume that one of the DWs is not moving 
(e.g., due to pinning). Let us assume that the classical par- 
ticle representing the other DW is located in one of the 
potential minima. If the distance between the minima is 
small, the particle can tunnel through the barrier, so that 
the other DW can be presented as a delocalized quantum 
particle. The Hamiltonian which describes it, has the form 

H = t X^i( c L f7 CTCT' c i+i-o'' +h.c.) where spin up and down 
states correspond to the collinear and anti-collinear orien- 
tations of the second DW with respect to the first one. 

Classically, if we consider a series of DWs pinned at 
a specific distance from each other, R*, inside a wire then 
the total energy of the system is given by an XY-model 

E = -J2\Ji cos(0 i+ i - 9i) + J 2 cos(0 i+2 - 9i)] , (11) 

i 

where only the nearest neighbour and next nearest neigh- 
bour interactions are included. Then we have J\ = 
-E int (0,R*) and J 2 = -E int (0,2R*). In general we 
can have J\ and J 2 as either negative or positive and ei- 
ther | Ji| > | J2I or IJ2I > I J\\. There are several possible 
set-ups admitting a simple solution. If |Ji| 3> IJ2I then 
we have either an antiferromagnet: 6>j + i — Oi = ir (for 




t - 1 th DW 



Figure 3 A schematic of the spiral structure of the DWs 
for J 2 < and \2J 2 \ > Ji|, the configuration of three of 
the DWs is shown. 9 = 9 i+1 - 6>; = cos" 1 [- Ji/2 J 2 }. The 
DW orientation, defined at the centre of the DW, and the 
bulk magnetization are shown. 

J 1 < 0), or a ferromagnetic arrangement (9; + i —6*^ = 
(for J 1 > 0). If we take Ji — > then we get two sublat- 
tices with either AFM or FM arrangements depending on 
the sign of J 2 . If we have 2J2 < J\ and 2J2 < — Ji, i-e. 
J 2 < and I2J2I > |Ji|, then the model is minimized 
by cos[6*.; + i — 9i] = — J\j2J 2 . In this case we obtain a 
spiral structure of the DW orientations through the wire, 
see figurte[3] Experimentally, similar spiral structures have 
been already observed I15II16I . 

3 Conclusion In summary we have found that there 
exists an RKKY-like electron mediated interaction be- 
tween DWs in a ferromagnetic sample that is long-range 
and oscillating. This interaction remains, independent of 
any current flowing through the system, in addition to pre- 
viously found results for current mediated interactions in 
nanowires. 
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